Cosmic ray electrons and positrons from supernova explosions of massive stars 
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We attribute the recently discovered cosmic ray electron and cosmic ray positron excess compo- 
nents and their cutoffs to the acceleration in the supernova shock in the polar cap of exploding Wolf 
Rayet and Red Super Giant stars. Considering a spherical surface at some radius around such a 
star, the magnetic field is radial in the polar cap as opposed to most of 47r (the full solid angle), 
where the magnetic field is nearly tangential. This difference yields a flatter spectrum, and also 
an enhanced positron injection for the cosmic rays accelerated in the polar cap. This reasoning 
naturally explains the observations. Precise spectral measurements will be the test, as this predicts 
a simple E~ 2 spectrum for the new components in the source, steepened to E~ A in observations 
with an E~ A cutoff. 

PACS numbers: 98.70.Sa, 96.50.sb, 97.20.Pm, 97.60. Bw 



Recently an excess of both cosmic ray (CR) positrons 
and cosmic ray electrons has been detected by three in- 
struments, PAMELA 0, ATIC 0, andH.E.S.S. 0. The 
ATIC and H.E.S.S. results on cosmic ray electrons are 
consistent with a discovery of an excess, compared to the 
normal measured spectrum of ^- 3 - 26±0 - 06 [j; W e em- 
phasize that also the H.E.S.S. result is above such an 
extrapolation. Both excesses take the form of a flatter 
component emerging from below a steeper, perhaps nor- 
mal, component; this has been interpreted in many ways, 
such as, e.g., the decay of a new particle representing dark 
matter [f|, or as evidence of a nearby special source Q. 
Here we wish to point out that such spectral components 
are expected from particle acceleration in the explosion of 
stars with magnetic winds 0] : There is a small polar cap 
component, where the acceleration in a supernova (SN) 
shock proceeds with the magnetic field parallel to the 
shock normal, yielding an E~ 2 spectrum [8j . At the same 
time, over most of Ait the magnetic field is best approx- 
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imated by an Archimedian spiral [9(, impl ying a near- 
perpendicular configuration for acceleration lOj . In such 
a situation the curvature becomes important, and the ac- 
celeration gives a slightly steeper spectrum and is faster 
- once injection has been effected. This was discussed 
again recently by Meli & Biermann [Til . The spectrum 
predicted is £-7/3-o.02±o.02 [^[^[ij The normal ob- 
served cosmic ray electron spectrum for E > 10 GeV 
provides a test, since it gives E~ 2 26±0 06 after correcting 
for losses [J, [l5[ . 

Our approach here is to adopt the following point of 
view [H Q : 

a) Most of the interactions of cosmic rays happen near 
the sources, and the escape from the Galaxy is governed 
by a simple Kolmogorov description [ljj • 

b) The bend in the spectrum, the knee at 10 15 eV is 
due to spatial limitations given by a shocked shell racing 
through a Parker- type magnetic wind [l2l . [l4| . 

c) We do invoke the physics of the stars that explode [l8| , 
and distinguish three zero age mass ranges of massive 
stars which explode: The stars between about 8 and 15 
solar masses, which explode into the interstellar medium; 
the stars between 15 and about 25 solar masses which 
explode as Red Super Giant (RSG) stars, and the stars 
above about 25 solar masses, which explode as Blue Su- 
per Giant, or Wolf Rayet (WR) stars. Both RSG stars 
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and WR stars explode into their stellar wind [l3|, Il9j ]. 
which is magnetic, and enriched from exposing the deeper 
layers of the star through mass ejections [20]. 

It was shown in Ref. [2l[ that such a combination of 
two components, a dominant E~ 7 ^ spectrum, with an 
additional E~ 2 , injected at the level of a few percent, 
yields a good fit to the cosmic ray air shower data. Pro- 
tons and heavier nuclei spectra are well described right 
through the knee of the cosmic ray spectrum (near 10 15 
eV) after accounting for transport that steepens the spec- 
tra by i? -1 / 3 . At this energy the polar cap component 
increases the overall flux by about a factor of 2, before 
cutting off due to spatial constraints [HI, H2] • The theo- 
retical assumptions about acceleration in highly oblique 
shocks 0, El used in RefJH were justified by the 
good agreement with data [23|: prediction E~ 8 / 3 , data 
^-2.68±o.o2 rpkj g ear jj er success encouraged us to apply 

precisely the same concept here. Additional enhance- 
ments of magnetic fields may occur also in such a situa- 
tion 01. 

We find, that using just the parameters of earlier pa- 
pers it is possible to explain both the cosmic ray electron 
spectrum and the cosmic ray positron excess. 

At GeV energies most of the cosmic ray (CR) elec- 
trons are accelerated in supernova (SN) shocks, running 
through the interstellar medium (ISM), the ISM-SN CRs. 
This predicts a spectrum of £- 2 -42±o.o4 [25], in agree- 
ment with radio data of other galaxies, for which the 
leakage energy dependence modifies the predicted spec- 
trum to £-2-75±o.o4 pe]. Already the data beyond about 
10 GeV suggest that the wind-supernovae cosmic rays 
may have taken over also for cosmic ray electrons [3] ■ 

Helium, Carbon and heavier nuclei give an indica- 
tion [lH of what fraction of 4-zr the polar cap component 
may have. This component reaches the same flux by itself 
as the rest of 4tt at 3 10 6 GeV/nucleus for CNO, yielding 
a surface fraction of about 2 percent. Since this fraction 
does not depend on distance from the star (see below) 
we assume that it defines the energy when the polar cap 
spectrum exceeds the E~ 7 1 3 component. 

Electrons, however, are injected at about 30 McV, the 
lowest energy at which they see the shock [23]. This 
number derives from the injection condition for electrons, 
that they must "see" the waves excited by the ions freshly 
injected by shocks in the assumption, that the plasma is 
dominated by ionized Hydrogen, and that the shock ve- 
locity is about 10,000 km/s. In a Wolf-Rayet star wind 
the main element is, however, not Hydrogen, but heavier 
nuclei, and already before the star explodes as a super- 
nova, there are accelerated electrons: The velocity of the 
shocks caused by instabilities in the radiation driving is 
of order > 1000 km/s [28[, and so the electron energy at 
injection then is at about > 6 MeV. 

This immediately implies that the polar cap compo- 
nent of cosmic ray electrons should rise to a flux equal of 
the rest at > 400 GeV, matching the uncertain observed 
energy of about 300 - 500 GeV. So at this energy the 
sum of the two components is twice the base spectral 



component. We interpret the ATIC data here as a £ 3 
component, rising above the base spectral component of 
£-10/3 arounc i 30 to 100 GeV. 

Cosmic ray positrons derive from collisions of nuclei, 
and formation of nuclei to the left of the valley of stability 
which decay in /3 + -emission, and also from pion produc- 
tion and decay. However, here we have to remember, 
that acceleration is faster for perpendicular shocks, by a 
factor up to c/(3V s h), probably more like 2-3 11]. This 
implies that the polar cap component is more efficient 
in producing positrons because of its slower acceleration 
and higher interaction probability. Since the hadronic in- 
teraction cross section is almost energy independent and 
does not introduce a break in the spectrum, the polar cap 
component becomes dominant at an energy between 2 to 
10 3 lower than for electrons, i.e. between 0.5 to 60 GeV. 
30 GeV seems to be compatible with the data, suggest- 
ing that the enhancement given by perpendicular shock 
acceleration is about a factor of 2 - 3. However, as there 
is a second source of positrons at lower energy, resulting 
from interaction in the immediate environment of mas- 
sive exploding stars and in the interstellar space [2^ . [30j | , 
the cross-over may be at lower energy, suggesting a pos- 
sibly higher efficiency enhancement. In a CR-positron to 
CR-(electron+positron) ratio this results in a rise with 

£+l/3_ 

From all these interaction sites here should be a cor- 
responding neutrino-emission with a spectrum of E~ 2 . 
On the other hand, as the cosmic ray electrons approach 
a spectrum of E~ 3 themselves, in the ratio positrons to 
electrons we approach a constant from somewhere in the 
range 30 - 100 GeV, when both electron and positron 
components are dominated by the polar cap. 

Here we discuss the second positron component at low 
energy, introduced above, which distorts the positron 
spectrum: 

The wave-field in the magnetic field excited by the cos- 
mic rays of spectrum E~ 7 ^ 3 [l3j] in the predecessor stellar 
wind naturally yields a specific spectrum of turbulence, 
in energy per volume and wave number I(k) ~ fc -13 / 9 , 
where k = 27r/r 9 , and r g = pc/ (ZeB), the Larmor 
radius (here p is the momentum of the particle, Z its 
charge, c the speed of light, e the elementary charge, and 
B the ambient magnetic field component perpendicular 
to the motion of the particle). This spectrum of mag- 
netic irregularities then governs the transport of cosmic 
rays and the cosmic ray interaction as a function of en- 
ergy. This in turn gives rise to a secondary to primary 
ratio going as E~ 5 / 9 (26|. This prediction was confirmed 
in Ref. [3l| which showed that the best fit of the second- 
ary to primary ratio had an energy dependence of _E~°' 54 . 
However, there is also another spectral component of tur- 
bulence induced by instabilities leading to many weak 
shock waves 7 (A;) ~ fc~ 2 , which is steeper. The total 
summed spectrum of turbulence has then a cross-over to- 
wards lower wave-numbers, corresponding to higher par- 
ticle energies; this spectrum I{k) ~ k~ 2 induces no en- 
ergy dependence of the production of secondaries. For 
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the most massive stars, exploding as WR stars, we esti- 
mate the cross-over to correspond to somewhere near 10 
GeV in electron energy. For somewhat lower mass stars, 
those exploding as red super giant stars, we argue, that 
there the cross-over between the two spectral regimes of 
turbulence is at lower energies, or higher wavenumbers, 
since the winds are less powerful. The cosmic ray in- 
duced turbulence is driven by the mass flow through the 
supernova shock, and so a wind of lesser density produces 
a weaker cosmic ray induced wave field. This results in 
secondaries having the same spectrum as the primaries, 
for WR stars above about 10 GeV, and for RSG stars at 
much lower energy, disregarding for a moment the polar 
cap component. 

This helps understand the diffuse gamma-ray emission 
of the disk of our Galaxy [32], as interaction near the 
RSG stars, more abundant than the WR stars. This 
model interprets all the secondary to primary ratios at 
low energy [31j . 

This may then explain the low energy cosmic ray 
positrons, seen with PAMELA [l[, as argued earlier. 

The acceleration time for cosmic ray particles in strong 
shocks is [1, E(| T acc — 8n/Vg h . Comparison of r acc to 
the synchrotron loss time T syn = 6irm e c/aTjeB 2 gives a 
limit on the maximal energy of the electrons or positrons. 

In the polar cap B ~ r~ 2 , while over most of Air B ~ 
r _1 . In the polar cap we assume maximal turbulence and 
Bohm diffusion k = (1/3)vlc. In the rest of the surface 
we adopt the approximation [ici| of n = V s h ■ The 
leads to maximum y e ,max = 1.5 • 10 6 BqI/^ V s h,g (r/r ) 
for the polar cap case, where B g is the magnetic field 
at radius 10 14 cm in units of 3 Gaufi, and V s h,a is the SN 
shock velocity in units of 10 9 cm/s. For most of 47r the 

maximum is j e ,max = 4.5 • 10 6 B^l'^ V^J (r/r ) 1/2 . 
Using the adopted values of the magnetic field and shock 
velocity, at 10 16 cm radial distance covered by the SN- 
shock in the wind the maximal energy of the CR electrons 
from the polar cap will dominate. Since the expressions 
used for k for both segments of the star's surface are 
simplified both energies are likely to be smaller. But 
Imax from the polar cap increases as r while for most of 
47r it increases as yfr [20|. This implies that the polar 
cap component will quickly pass the other component in 
maximal energy. 

Last, we check whether the angular fraction is really 
distance independent, as required by such a model. Com- 
bining the results derived in Refs. [HI EH, EH EH EH we 
match the acceleration time scale in the polar cap with 
the acceleration time over most of 47r. In the limit of a 
small angular extent of the polar cap this gives: 8 = 
3 (B ,p. c ./B 0A7r ) (V sh /c) suggesting for V sh /c ~ 0.03 a 
ratio of surface magnetic fields of Bq iP . c J Bqajt — 0.6, 
very close to unity and independent of radius. It is an 
interesting question if for some stars this ratio might be 
different. Accelerated electrons have energy loss given 
by the synchrotron loss time T syn and they can also es- 
cape the Galaxy with r/ ea fc ~ 7 e ~ 1 / 3 . Then we have 



N(E) ~ E 7 / 3 1 spectrum in the limit T syn < Ti ea k 
and N(E) ~ E^ 7 / 3 ^ 1 / 3 , in the opposite limit. We can 
estimate from observations, radio measurements as well 
as direct data, that the switch-over is near 10 - 20 GeV. 

CR electrons reach us in a random walk, in which 
the distance actually travelled is given by r ~ Ar \/~N, 
where Ar(E) is the diffusion scattering mean free path, 
the step, and N is the number of uncorrelated steps. 
The time this takes is given by t ~ (Ar/c)N, while 
Ar ~ E 1 / 3 and the timescale t is ~ 1/E. There- 
fore r ~ E- 1 / 3 . The region from which we obtain CR 
electrons is the volume for which this time is less than 
the synchrotron time T syn . The source volume is propor- 
tional to r 3 ^ E^ 1 . Adopting the viewpoint that at high 
energies the polar cap component of the CR contribution 
of the the wind-SNe dominates implies then a combined 
spectrum of ~ E~ 3 ~ x . So the CR electron spectrum is 
predicted in this approximation to be E~ 10 / 3 , then E~ 3 , 
and thereafter E~ 4 which is consistent with the H.E.S.S. 
data Q. Nearby massive stars exist [33| and the associ- 
ated supernovae [34| may have provided a large fraction 
of the cosmic rays we observe. A test for our predictions 
would be a measurement of the exact spectrum of both 
cosmic ray electrons as well as cosmic ray positrons. This 
model simply predicts that their spectrum is a simple ad- 
ditional E~ 2 component at source. 

The positron fraction should reach a plateau, the ex- 
act number depending on the dominant path to pro- 
duce positrons, decay from isotopes pushed by photo- 
dissociation and spallation off the valley of stability in a 
N, Z plot, or just simply pion decay. 

Ion collisions produce very few if any anti-protons and 
contributions from WR star explosions will be negligi- 
ble. The RSG explosions will produce anti-protons [35j . 
and their polar cap contribution should come up at some 
higher energy, perhaps above 30 - 100 GeV. 

Obviously, both cosmic ray electrons as well cosmic ray 
positrons are in their respective loss limit [l5| at such en- 
ergies, above about 10 - 20 GeV, where losses overpower 
diffusion and steepen the spectrum by unity; the limited 
spatial reach gives a further steepening by E~ l . There- 
fore there must be some yet higher energy where the 
losses from the nearest most recent source cut everything 
completely off; the ATIC and H.E.S.S. data suggest that 
this happens beyond energies of several TeV. It would be 
very interesting to measure the positron component to 
this energy. 

We have proposed a simple explanation for the cos- 
mic ray electron and cosmic ray positron components, 
in terms of the magnetic field topology in a magnetic 
wind 0, and SN- induced shock acceleration in such a 
topology [lflEHl- The new component is just that popu- 
lation of energetic particles accelerated in the polar cap of 
massive magnetic stars with winds, when they explode. 
We attempted to explain a) the low energy PAMELA 
data, b) the higher energy PAMELA data, c) the ATIC 
data, and d) the H.E.S.S. data, all in the context of a 
basic picture proposed and worked out earlier [IH [H|, 



4 



+ 





... M, 




dN/dE^ E~ 


0/3 ^--3 i Iv 




• H.E.S.S. (2008), 


vith 1 5 % correction 


\ 1 


■ ATIC (2008) - po 


n:s excluded from modeling 




■ ATIC (2008) 

1 r , ., , I , , , , 1 , , , r 1 , 


, , i , , , , i , , , , i , , , , i , , 




0.5 0.5 1 


1.5 2 2.5 3 


3.5 




E -/3 










1 , , \ , 


T Pamela 




(2008) 



-0.5 0.5 V 7.5 2 2.5 3 3.5 4 



log(E/GeV) 

FIG. 1: Model prediction with ATIC and H.E.S.S data for 
cosmic ray electrons: E~ 10 / 3 , E~ 3 , and E~ 4 (upper panel), 
and for PAMELA data on the positron/electron ratio: E~ 5 ^ 
and E 1 / 3 (lower panel). See Ref. [3(| for charge dependent 
solar modulation arguments. 

and consistent with other measurements. 

In the cosmic rays produced by a shock running 



through such a wind there is always a polar cap com- 
ponent, of a few percent strength at injection, but with 
a flatter spectrum. 

The explosions of Wolf Rayet stars and their cousins, 
the Red Super Giant stars, into their respective magnetic 
wind, may play a key role in allowing us to understand 
the physics of cosmic rays. 
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